In-plane anisotropy of the electronic structure for the charge/orbital-ordered state in 

half-doped manganite with layered structure 
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We report on the in-plane anisotropy of the electronic response in the spin/charge/orbital or- 
dered phase of a half-doped layered-structure manganite. The optical conductivity spectra for a 
single domain of Eui/2Ca3/2Mn04 unambiguously show the anisotropic charge dynamics which well 
corresponds to the theoretical calculation: the optical conductivity with the polarization along the 
zigzag ferromagnetic chain direction exhibits a smaller gap and a larger intensity at lower ener- 
gies than that of the perpendicular polarization mostly due to the charge/orbital ordering and the 
associated quantum interference effect. 

PACS numbers: 75.47.Gk, 75.30.Et, 78.40.Ha 



Recent studies on strongly correlated electronic sys- 
tems have revealed many unexpected phenomena due to 
interplay among the spin, charge and orbital degrees of 
freedom. In particular, nano-scale self-organization and 
the associated complexity of electronic systems are issues 
of current interest. One of the representative phenomena 
is the spin/charge/orbital ordering in transition metal ox- 
ides (TMOs), which governs the electrical, optical, and 
mechanical properties of the system [lj . Doping carriers 
into the parent Mott insulating compound causes a va- 
riety of ordering phenomena including one-dimensional 
(ID) stripe formation as observed in some layered Cu- 
and Ni-oxides. In Mn-oxides near half-doping there ap- 
pears a cooperatively ordered pattern of charge, spin, and 
orbital, referred to as CE-type magnetic order coupled to 
the charge/orbital ordering (CO-OO) HHSH0, which 
is involved in the metal-insulator transition with mag- 
netic field and other various control parameters Q. Be- 
low the CO-00 transition temperature Tco-oOi charges 
are ordered in a checkerboard-type with ideally alternat- 
ing Mn 3+ /Mn 4+ sites, i.e., the 3d e 4 /e° configuration 
in addition to localized tig spins, and e g orbitals are 
aligned in the zigzag chain type. The direction of e g 
orbitals defines the zigzag chain and stripe [Fig. 1(c)]. 
Below the Neel temperature T/v, which is usually far be- 
low Tco-oo, the spins of 3d electrons at Mn sites are 
ordered ferromagnetically along the orbital zigzag chain 
path, while antiferromagnetically between the chains. 
This well-organized charge/spin/orbital ordering occurs 
in the ab plane for a variety of manganites including lay- 
ered compounds, and generically has the effective two- 
dimensional nature in electronic response even for the 
pseudo-cubic perovskite 

To access the attribute of the charge/orbital order- 
ing in more depth, we pay attention to the in-plane 
anisotropy of the electronic structure in the CO-00 
phase, in particular, the optical response. The spin order 



may invoke the simplest trivial source of the anisotropy. 
Namely, well below TV where the ferromagnetic chains 
are coupled antiferromagnetically, the local electron hop- 
ping between the adjacent chains costs a larger energy 
typically of the order of JhS (JptHund's rule coupling, 
S: magnitude of the ti g spin), than that along the zigzag 
chain path. Since the optical processes occur locally, they 
are insensitive to the range of the magnetic ordering, and 
the anisotropy remains even above Tjy. This scenario, 
however, is oversimplified and may not apply directly to 
real experimental systems that show the drastic change 
of electronic response below Tco-oo^ which is even well 
above T/v- Therefore, we must seriously consider the ef- 
fects of CO-00 on the anisotropy of the low-energy op- 
tical response. 

To elucidate the effects of CO-00 in the in-plane op- 
tical response, we take the following minimal model for 
a ferromagnetic zigzag chain in the CE phase; 

r,r' a,/3=l,2 r 

+Ej T X (ul + vl - 2u r (n r i - n r2 ) - 2u r (4 1 c r2 + 42 c riX] 
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With the annihilation (creation) operators c ra (cj a ) and 



the number operator n ra = cj 



of the electron at site 



r with the e g orbital a (—1 or 2 for d 3z 2_ r 2 or d x 2_ y 2 
respectively), the Hamiltonian given by Eq. Q contains 
the orbital-dependent electron transfer t^f between the 
adjacent sites r and r' with the same spin polarization 
as the first term, and the local Coulomb and the Jahn- 
Teller interactions U and Ejt as the second and the 
third terms, respectively. u r and v r denote the local 
Jahn- Teller modes for the MnOg octahedron centered at 
the Mn site r, as shown in the inset of Fig. 1 (b). We 
have neglected the phonon dynamics of these Jahn- Teller 
modes. This frozen Jahn- Teller distortion corresponds 
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FIG. 1: (Color online) (a) Optical conductivity spectra along 
the chain and the stripe directions from (a) the calculation us- 
ing the Hartree-Fock method and (b) the experiment for the 
single domain of ECMO. In the inset of (a) the result of the 
exact diagonalization study is displayed. In the inset of (b) 
two Jahn- Teller modes, u and v, are drawn with a MnOe octa- 
hedron, (c) Schematic diagram for the charge/orbital ordered 
state at half-doping. Right bottom panels depict difference in 
the phases for the local hopping and current matrix elements 
from the d 3x 2_ r 2 /d 3y 2_ r 2 orbitals at the electron-rich sites to 
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(d) the d x 2_ y 2 and (e) the d iz 2_ 
poor ones. 



orbitals at the electron- 



to the Frank-Condon approximation which is justified 
for the strong electron-phonon coupling case as in the 
present situation. As for the orbital and directional de- 
pendencies of the transfer integrals, we take the Slater- 
Koster's value with Vdda — t — 0.25 eV, which is taken 
as the energy unit. 

We solve this Hamiltonian within the Hartree-Fock ap- 
proximation to obtain the self-consistent values for the 
expectation values {c\ a c r ' a i) as well as the Jahn- Teller 
distortions u r and v T . We note that Hartree-Fock treat- 
ment provides a reasonable description of the ground 
state even at the large U. Since the result is insensitive to 
the value of U, we may take a typical value for transition- 
metal ion, U = 8t. In the present case of half doping, we 
have obtained the CO-00 pattern characterized by the 
charge disproportionation Sn and the dominant orbital 
character of d 3x 2 _ r 2 / d 3y 2 _ r 2 at the electron-r ich (i^a) 
sites and d x 2_ 2 at the electron-poor ( 1 ~ (5n ) ones. For 



this ordering the Jahn- Teller coupling is indispensable, as 
previously pointed out |10| , and the ground-state proper- 
ties are insensitive to the choice of U. Then, we employ 
the value Ejt = 2t to reproduce the experimental result 
of the optical gap. This choice of Ejt = 2t gives the 
charge disproportionation of Sn = 0.56, which is inter- 
mediate between the uniform value Sn = for Ejt = 
and the ideal value Sn = 1 for Ejt —* 00. 

Using the Hartree-Fock solution, we have calculated 
the optical conductivity <j cham and <j stripc along the chain 
and the stripe directions, respectively, with the corre- 



sponding current operators J cham 
jstripc = ( jx _ w here 



(J x + J y )/V2 and 
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J x ' v = — -7 



N ^ dk x 

k,a,/3 



K' t 

c k,a C k,/3 



(2) 



with t 

_cha 



a,/3 



E«Cr+5 r ' k ' S - The calculated results for 
and <r stripe in the case of Ejt = 2i are shown 
in Fig. 1(a). This calculation reveals that the optical 
conductivity along the chain direction exhibits a smaller 
optical gap and a larger spectral weight distribution at 
lower energies than along the stripe direction. We have 
also done the exact-diagonalization study of the same 
model with the Jahn- Teller distortion being fixed at the 
ground state value (Frank-Condon approximation), and 
found that the global structure of cr cham and cr stripe and 
the emergence of the anisotropy do not alter even be- 
yond the Hartree-Fock, as shown in the inset of Fig 1(a). 
The anisotropy originates from the difference in the phase 
factors for the combination of J x and J y in the orbital- 
ordered state and is an increasing function of the CO- 
OO order parameters. More explicitly, the local hopping 
process from the d^ x 2_ r 2 /d 3y 2_ r 2 orbitals at the electron- 
rich sites to the d x 2_ y 2 orbital at the electron-poor site, 
which is shown in Fig. 1(d), exhibits additive and neg- 
ative interference in the chain and the stripe directions, 
respectively, leading to the strong anisotropy, while the 
hopping to the d 3z 2_ r 2 orbital shown in Fig. 1(e) does 
not. 

To address the in-plane anisotropy experimentally, the 
measurement for an in-plane twin-free single domain 
sample is highly required, but has not been performed 
so far due to a great difficulty to obtain a sample large 
enough for optical measurement (~ 100 fim in case of 
using microscope probe). To overcome this technical 
problem we chose a 2D single layered structure which 
is quite suitable for probing the Mn-0 plane where the 
CO-00 occurs. Due to the two-dimensional electronic 
nature of the charge/orbital ordering concerned here, our 
study with layered manganites is directly applicable to 
the case of the pseudo-cubic perovskite. For the well- 
studied half-doped material Lai/ 2 Sr 3 / 2 Mn04 (LSMO), 
a typical domain size was found to be ~ 10 x 10 fim 2 
|ll| . not sufficiently large for our purpose. Consider- 
ing that the lattice distortion by the small A-site ions 
helps the domain size larger we newly synthesized the 
Eu 1 / 2 Ca 3 / 2 Mn04 (ECMO) compound where the aver- 
age size of A-site ions (~ 1.165 A) is much smaller than 
that in LSMO (~ 1.285 A) [l2j. The T CO -oo for the 
ECMO compound is as high as ~ 325 K perhaps due to 
the narrower bandwidth, which is identified from anoma- 
lies in resistivity and magnetization, and a superlattice 
peak in the electron diffraction pattern (EDP). The Tjv is 
estimated to be 120 - 130 K far below T c0 -00 E3- Us- 
ing the polarized microscope we could identify the stripe- 
type domain pattern in the cleaved ab surface. While the 
typical widths of domains are 10-50 /im, we succeeded in 
finding an exceptionally large single domain with ~ 300 
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FIG. 2: Reflectivity spectra (left panels) and optical conduc- 
tivity spectra (right panels) of Eui/2Ca 3 / 2 Mn04 along the a' 
and b' axes. The panels at each temperature are shifted up in 
sequence. The 335 K data are included in the top panels to 
clarify the temperature-dependence. The dotted line in the 
right panels is for guidance. 



fim. width, which enabled us to perform the polarization- 
dependent optical measurement for the CO-00 phase. 

A Euiy2Ca3/2Mn04 single crystal with half-doping was 
grown by the floating zone method. For the reflectivity 
measurement of a single domain, we used the microscope- 
equipped spectrometers with the beam size, ~150 x 100 
/im 2 . With proper polarizers the probed light is polarized 
along the optical axis (a' or &'), which is rotated by 45 
degree from the crystal axes of the tetragonal cell. The 
EDP experiment on this material has revealed that the a' 
(&') axis chosen for our measurement should correspond 
to the chain (stripe) direction. T-dependent reflectivity 
spectra R(u>) at nearly normal incidence were measured 
in a photon energy region from 0.06 to 5 eV with the tem- 
perature variation from 10 to 430 K. Since the anisotropic 
feature in R(lo) is terminated in the charge transfer exci- 
tation region above 2.5 eV, T-dependent R(lu) below 5 eV 
are connected smoothly to room temperature R{uj) above 
5 to 40 eV obtained with a twinned sample. Optical con- 
ductivity spectra c(w) were obtained from the measured 
R(uj) using the Kramers-Kronig transformation. 

Figure 2 shows the reflectivity and the optical conduc- 
tivity spectra with polarization dependence at various 
T. Both of optical spectra along a' and b' axes show a 
similar trend with T, exhibiting significant changes re- 
lated to the CO-OO. As T is lowered below Tco-oo, 
the reflectivity level is suppressed in a low energy re- 
gion below ~ 1 eV, while enhanced in the higher energy 
region. The cr(u>) calculated from the measured i?(w) 
exhibits that a distinct absorption below 2 eV, which is 
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FIG. 3: (a) Optical gap, 2A opt and (b) absorption peak below 
2 eV, Wpeak along the a' and b' axes, (c) Differential conduc- 
tivity Aa(u}) — a a i (u>) — a b i (lj) at various temperatures. Inset 
displays the temperature dependence of AN e ff(u) = 1.3 eV). 
The unit is f2 _1 cm -2 . 



associated with the electron hopping from Mn 3+ (occu- 
pied e g ) to Mn 4+ (unoccupied e g ) [14(, shifts to a higher 
energy with the sharper spectral shape at lower T be- 
low Tco-oo, leading to a larger optical gap. The overall 
spectral feature is in accord with the previous results for 
twinned LSMO samples ITD . Il5| . and the related theo- 
retical calculations [la. fl7| . Apart from this common 
T-dependence, we could discern a clear development of 
anisotropy in optical spectra below Tco-oo- The reflec- 
tivity spectra at 335 K (above Tco-oo) appear to be 
nearly identical for both polarizations. In contrast, at 
320 K (below Tco-oo) the distinct anisotropy is identi- 
fied; the level of reflectivity in the b' axis is lower below 
1.3 eV, while higher above 1.3 eV, than in the a'-axis. 
The difference between R(u>) along the two axes is larger 
at lower T. Similarly, the spectra polarized along 

the 6'-axis are distinguished from those along the a'-axis 
below Tco-OOj exhibiting the higher energy of absorp- 
tion band and the larger optical gap. 

To clarify the in-plane anisotropy, we display the po- 
larization dependence of the optical gap, 2A opt , and the 
absorption position below 2 eV, cj pca k, in Figs. 3(a) and 
3(b). The 2A opt values were estimated from crossing 
points of abscissa with linear extrapolations of (t(oj). At 
T < Tco-oo both the quantities have larger values along 
the 6'-axis than the a'-axis. The anisotropic difference 
in 2A opt and a> pe ak are estimated to be ~ 0.12 eV and 
~ 0.35 eV, respectively, at the lowest temperature, and 
hence the degree of anisotropy is as large as 20 - 30 %. 
It is noted that such sizable anisotropy cannot originate 
from the orthorhombic crystal structure alone. While the 
phonon spectra reveal the lowering of crystal symmetry 
below Tco-oo 0, 0> the degree of orthorhombicity 
is too small (b'/a' — 1 < 0.01) to give rise to the sig- 
nificant electronic anisotropy observed here. In fact, the 
observed anisotropy above Tco-oo-, which is supposed 
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to come from the inherent lattice orthorhombicity apart 
from CO-OO, is extremely small. 

Notably, the electronic anisotropy unveiled by the ex- 
periments is in good agreement with the theoretical pre- 
diction as shown in Fig. 1. Because the a' (b 1 ) axis cor- 
responds to the chain (stripe) direction, our combined 
experimental and theoretical studies unambiguously re- 
veal the formation of the smaller energy gap state along 
the chain direction. This finding could be a strong evi- 
dence of more effective electron hopping along the chain 
direction in the CO-00 phase. Clearly the in-plane 
anisotropy identified below Tqo-OO is attributed to the 
unidirectional zigzag chain arrangement of the e g orbitals 
accompanied by the charge ordering. The orbital order- 
ing and the associated anisotropic transfer integrals lead- 
ing to the interference of various hopping processes play 
a crucial role in the electronic structure of the CO-00 
phase. The direction of orbital stripe (or orbital chain) 
is closely related to that of the optical anisotropy. The 
different optical gap value revealed in our spectra is ex- 
pected to lead to the high anisotropy in dc conductiv- 
ity. Applying an external electric field might change the 
direction of the orbital stripe and the corresponding di- 
rection of optical anisotropy. In this sense the orbital 
stripe might work as a parameter of a sort of optical 
polarization-related device. 

Finally we discuss the correlation between the 
anisotropy and the magnetic correlation by tracing their 
T-dependent evolutions. To get some insight into the 
anisotropic feature in <t(uj) below 2.5 eV, we estimated 
the effective spectral weight by integration of a(u>), 
N e g(uj c ) — J Q C duj'a(uj'), and found that N g(u] c = 2.5 
eV) are nearly identical in both axes. This indicates that 
the optical anisotropy observed in the CO-00 state is as- 
sociated with the distribution of the local spectral weight 



below 2.5 eV. To clarify this, we display the differential 
conductivity spectra Act(cj) = a a i(uj) — crj'(w) at vari- 
ous T in Fig. 3(c). While the Aer(w) above Tco-oo 
exhibits quite a small spectral weight with negligible In- 
dependence, the distinct two peaks with opposite signs, 
which originate from the accumulation of lower-energy 
spectral weight along the a 1 axis, develop significantly 
below Tco-oo- We also estimated the differential effec- 
tive spectral weight AN c g= f^ a duj' Act (a/) with u c = 1.3 
eV, where the ct(oj) in both axes coincide at the lowest 
T. As shown in the inset of Fig. 4(c), AN e g exhibits a 
sizable value with a steep increase below Tco-oOi an d 
then appears to be rather saturated at lower T, not ex- 
hibiting a distinct change near Tn- From this it appears 
that the three-dimensional (3D) long-range spin order- 
ing does not affect the optical anisotropy significantly. 
When the orbital ordering occurs, the short-range mag- 
netic correlation could be presumably formed along the 
zigzag chain even above TV, which have been identified 
by the neutron diffraction measurement for some CO- 
OO manganites El- No significant change of the 
anisotropy near T/v implies that this ferromagnetic fluc- 
tuation along the chain could be effective for the forma- 
tion of the quasi-lD electron hopping, irrespective of the 
onset of 3D long-range magnetic order. This may be in- 
dicative of the predominant role of the orbital ordering 
in the magnetic correlation. 

In summary, we have presented the experimental 
and theoretical studies on the optical spectra in the 
charge /orbital ordered plane with polarization depen- 
dence. The in-plane anisotropy is attributed to the pe- 
culiar orbital arrangement accompanied by the magnetic 
correlation. 
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